Suppression of the "quasiclassical" proximity gap in correlated-metal— superconductor 

structures 
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We study the energy and spatial dependence of the local density of states in a superconductor- 
correlated-metal-superconductor Josephson junction, where the correlated metal is a non-Fermi 
liquid (described by the Falicov-Kimball model) . Many-body correlations are treated with dynamical 
mean-field theory, extended to inhomogeneous systems. While quasiclassical theories predict a 
minigap in the spectrum of a disordered Fermi liquid which is proximity-coupled within a mesoscopic 
junction, we find that increasing electron correlations destroy any minigap that might be opened in 
the absence of many-body correlations. 
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Fermi-liquid metals have excitation spectra that typi- 
cally vary on energy scales of electron volts. Metals that 
become superconducting, have all low-energy electrons 
form pairs. Since it takes an energy equal to the super- 
conducting gap A to break a pair of electrons and form 
an excitation, there is a low-energy gap (on the order 
of meV) in the single-particle density of states (DOS) 
0. The original states at the Fermi level are "pushed" 
to excitations near ±A, which yields a singularity in 
the DOS at zero temperature (and a large peak at fi- 
nite T). When a superconductor (S) is connected to a 
normal metal (N) to form a SNS Josephson junction, 
the superconductivity leaks into the normal metal via 
the proximity effect and a weak link is established 
between the two S through the N. What happens to 
the low-energy electrons in the N is quite interesting. 
An electron near the Fermi level of the N is localized 
within the N because there are no single-particle states 
at low energy for it to scatter into within the S. In- 
stead, the electron is retroreflected into a hole in the N, 
and creates a superconducting pair in the S via a pro- 
cess called Andreev reflection This reflection occurs 
at the SN interface and the NS interface, creating An- 
dreev bound states with well-defined energy levels. These 
states are doubly degenerate, one carrying supercurrent 
to the right and one to the left. The states are broadened 
into peaks in the DOS when one averages over all dif- 
ferent perpendicular momenta. Thus, Andreev reflection 
mixes the electron and hole states in the same proportion 
that they are mixed to form Bogoliubov quasiparticles in 
the S, with weights determined by the self-consistency 
condition. Such partially superconducting properties of 
the proximity-coupled normal metal are responsible for 
Josephson effect, as well as other peculiar phenomena 
in inhomogeneous systems which have been drawing in- 
creased attention over the past decade due to advances 
in mesoscopic superconductivity 

Another aspect of the proximity effect is the modi- 
fication of the DOS in both the S ("inverse proximity 



effect " nil) and the N side of a SN boundary, which 
becomes most conspicuous in mesoscopic confined ge- 
ometries 0. For example, in a sufficiently long SNS 
junction at low enough temperature, the proximity of 
the superconductor induces a minigap in the local den- 
sity of states (LDOS) inside N interlayer that 
has chaotic classical dynamics [| 10 . The minigap is 
of the order of the Thouless energy i?Th — ^/idwcii, 
where tdwcU is the typical time spent by an electron dur- 
ing its diffusive motion |]7|,p|,pl| (idwoii — L'^/V for a N 
strongly coupled to a S, where V is the diffusion con- 
stant) or during its chaotic ballistic motion [|lO| through 
the N region of size L, before escaping into the super- 
conductor (for integrable classical dynamics in the TV, 
the LDOS is nonzero, but vanishes nearly linearly at 
the Fermi level |jl^). These results were obtained us- 
ing either quasiclassical approaches ||^-p|, pT[ o r mean- 
field treatment by random matrix theoryip^ . How- 
ever, recent calculations which include additional 
quantum effects through a supersymmetric non-linear a- 
model (NLSM) 0, show that mesoscopic fluctuations 
cause an exponentially small smearing of the quasiclas- 
sical gap in a diffusive SNS junction (the DOS tails 
appearing below the quasiclassical gap edge are due to 
prelocalized states jl4| which couple weakly to the S 
leads). Using tunneling spectroscopy, the proximity af- 
fected LDOS can be measured as a function of the dis- 
tance from the SN interface in both the normal metal jl^] 
and superconductor ||]. 

While both quasiclassical and NLSM calculations rely 
on the picture of well-defined (noninteracting) quasipar- 
ticles, little attention has been paid to proximity effects 
in systems where the N layers are dominated by strong 
many-body correlations (except in one dimension where 
anomalously enhanced DOS have been found in a Lut- 
tinger liquid coupled to a superconductor [l^). Here 
we explore the LDOS, in both the superconducting and 
correlated metal sides of a superconductor-correlated- 
metal-superconductor (SCmS) Josephson junction. The 
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FIG. 1. Local DOS on the normal plane closest to the 
center {z = 0) of SCmS Josephson junctions of thickness: 
(a) L — 5a and (b) L — 10a. The Cm interlayer is 
a non-Fermi liquid described by the FK model where the 
bulk DOS is flat in the plotted energy range around the 
band center ui — [e.g., N{uj) = 0.127 at Ufk = 2, see 
Fig. ^. The minigap in the LDOS of the Cm is open only 
for small enough i7pK, and does not scale as ^ which 
would be the quasiclassical prediction, but seems to be just 
a remnant of the minigap ~ A^//^ [inset in panel (a), where 
A(;7 = -2)/A(t/ = -1.7) f» 1.9] in a clean SNS junction 
which is gradually destroyed by increasing the strength of 
many-body correlations. 

Cm region is a non-Fermi liquid modeled by a Falicov- 
Kimball (FK) Hamiltonian [^. We find that increas- 
ing electron correlations completely destroy any initially 
open "quasiclassical" minigap in the DOS of the corre- 
lated metal, as shown in Fig. 0; this occurs due to the 
extensive broadening of the "Andreev bound states" by 
the scattering in the Cm. The appearance of a nonzero 
LDOS inside the S region within a distance on the or- 
der of the superconducting coherence length is plotted 
in Fig. ||. Our analysis is fully self-consistent (i.e., we 
take into account the suppression of the superconducting 
order parameter inside the S leads) at zero Josephson 
current across the junction. The calculation includes all 
(many-body) quantum effects encompassed in the 



FIG. 2. Local DOS in the superconducting side of the 
SCmS junction with thickness L — 5a and the Cm inter- 
layer described by the FK model with Ufk = 1.0. The planes 
are labeled by their distance from the center {z = 0) of the 
Cm region, in units of the superconducting coherence length 
~ 4a (the SCm interface occurs at z = 0.75^s). The re- 
gion inside the semi-infinite S leads where the self-consistent 
calculation is performed extends to z = 8.25^s. 

dynamical mean- field theory [ p^ , which has only re- 
cently been generalized to treat inhomogeneous normal 
systems and Josephson junctions | |20[ |. 

Thus, our principal result is substantially different 
from the standard lore of a proximity-induced "hard 
minigap" (i.e., no states inside an energy interval Eg 
ETh)i which is supported by both quasiclassical calcu- 
lations |^J^,|ll| and the picture of bound states induced 
by Andreev reflection [^. Such qualitative considera- 
tions [^ give an estimate for Eg, which comes close to 
the values of the quasiclassical minigap obtained from the 
solution of the Usadel equation: Eg = 0.78i?Th in the dif- 
fusive N layer of an INS structure ||T^ (/ is an insulator 
that specularly reflects), or Eg = 3.12i?Th in a diffu- 
sive SNS junction l^. The perturbative analysis of the 
electron-electron interaction, within quasiclassics, only 
generates a slightly smaller Eg The NLSM calcula- 
tions [|l3| for a SNS junction, with N being a disordered 
noninteracting electron system, flnd nonzero DOS at all 
E < Eg, but the tail of subgap states is small for good 
metals (except near Eg). Our results are complementary 
to these theories, showing how the minigap can disap- 
pear at a critical correlation strength, and are of special 
interest in understanding the limitations of a phenomeno- 
logical application of standard proximity-effect theory to 
experiments dealing with unconventional inhomogeneous 
structures, such as high- Tc Josephson junctions with un- 
derdoped cuprates (a strongly correlated electron sys- 
tem) playing the role of the "normal region" pl[ |. 

The SCmS Josephson junction is modeled by a Hamil- 
tonian 
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on an infinite set of stacked square lattice planes, whose 
connectivity is the same as a simple cubic lattice (with 
lattice constant a). Here (cio-) creates (destroys) an 
electron of spin a at site i, Uj — t (the energy unit) 
is the hopping integral between nearest neighbor sites i 
and j (both within the planes and between planes), Ui < 
is the attractive Hubbard interaction for sites within 
the superconducting planes, is the FK interaction 

for planes within the Cm region, and Wi is a classical 
variable that equals 1 if an A ion occupies site i and is 
zero if a i? ion occupies site i. The chemical potential 
ji is set equal to zero to yield half filling in the S and 
Cm. The negative-C/ Hubbard term describes the real- 
space pairing of electrons due to a local instantaneous 
attractive interaction. This generates a superconducting 
order in the S leads which, when treated in the Hartree- 
Fock approximation, is equivalent to conventional BCS 
theory, except that here the DOS is non-constant and 
provides the energy cutoff. 

The superconducting layers have Ui = —2 and Wi = 
for all sites. Such a homogeneous bulk superconduc- 
tor is characterized by the usual BCS parameters: the 
transition temperature Tc — O.llf, the zero-temperature 
order parameter A = 0.198t, and the coherence length 

= ftfF/(7rA) K, 4a. The Cm interlayer is described by 
a half-filled FK model in the symmetric limit of half fill- 
ing for the "ions" (wi) = 0.5. One can view the FK metal 
as a binary alloy of A and B ions at 50% concentration 
with Ufk being the difference in site energy between the 
A and B ionic sites. The many-body problem is solved by 
taking an annealed average that yields the coherent po- 
tential approximation. This is the simplest many-body 
problem which, nevertheless, mimics a metal-insulator 
transition (MIT) of the type seen in the repulsive Hub- 
bard model (except the metallic phase is not a Fermi 
liquid here) . In the bulk, the FK correlated metal under- 
goes a MIT at UpK ~ 4.9t (which is close to half of the 
bandwidth 6t). This is illustrated in Fig. || which shows 
the DOS in a bulk Cm as a function of Ufk- The DOS is 
independent of temperature ||2^ . Since the system is not 
a Fermi liquid for nonzero Ufk, the DOS first develops 
a pseudogap, and then is suppressed entirely to zero as 
the correlations increase, eventually driving the system 
into a correlated insulator. The opening of the gap is 
continuous. In order to focus only on the modification 
of the DOS induced by the proximity effect, we choose 
Ufk < 2 for the strength of Coulomb interaction in the 
FK correlated metal, which ensures that the bulk DOS 
around the band center is essentially constant. 

The problem of inhomogeneous superconductivity 
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FIG. 3. Electronic DOS (per spin) for the bulk Cm de- 
scribed by the FK model on a simple cubic lattice in the local 
approximation. The value of Ufk ranges from 1 to 7 in steps 
of 1. As Ufk increases, the DOS first develops a pseudogap 
and then a real gap. We use Ufk < 2 for the Cm region of 
the SCmS Josephson junction (Fig. |l]) where the bulk DOS 
is constant around half-filling. 

is solved by employing the Nambu-Gorkov matrix for- 
mulation for Green functions with a local self-energy . 
We treat the problem self-consistently in the complex 
order parameter for the part of the junction com- 
prised of the normal region (containing 5 or 10 planes) 
and 30 superconducting planes on each side of the Cm 
interlayer. Inside this "self-consistent-part" of the infi- 
nite SCmS junction, the superconducting gap Aj heals 
to its bulk value since all signatures of the inverse prox- 
imity effect are gone on the length scale of few ^5 away 
from the SCm interface. The calculation is performed at 
the temperature T = 0.09Tc where the BCS gap is fully 
developed. Details of our computational algorithm have 
been given elsewhere |2^]. The final result is the self- 
consistent Green function which allows us to compute 
the (many-body) LDOS as a function of plane position 
Zi and frequency lu from the real-axis analytic continua- 
tion, p{uj, Zi) = — ImG(ti; -I- iS; Zi, Zi)/!:. 

The strength of the superconductivity in the N inter- 
layer is quantified by the nonzero pair amplitude F{zi) = 
Ai/|[/i| (a two-particle property) which decays exponen- 
tially, due to the absence of an attractive interaction, on 
the length scale = ^/Wofh^k^ @ (or as a power 
law at zero temperature in a clean normal metal |p3| ). 
This is the length scale over which two thermal elec- 
trons with energy lo « TrfcsT in the TV, correlated by An- 
dreev reflection over the length scale L^j = ^/KD/2uj, lose 
their relative phase coherence (which then determines the 
thermodynamic critical Josephson current Ic). However, 
single-particle properties (like the DOS) ^,|l^, or kinetic 
quantities , can be influenced on a much longer length 
scale (where Josephson coupling vanishes) ^ S^jq for 
low-energy electrons lo <^ kT, which is ultimately limited 
by the mesoscopic phase-breaking length There- 
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fore, the anomalies in the DOS for low uj extend up to the 
energy dependent distances ^ from the SN in- 

terface (which is smeared upon approaching |[T5[|). In 
a "closed" geometry, where a finite-size N is disconnected 
from any electron reservoirs 0, this leads to a position 
independent minigap edge at Eg ~ £^Th for sufficiently 
long i<^Th ^ A diffusive interlayer [l <C L, where I is 
the mean free path). Thus, both the quasiclassical and 
NLSM descriptions of the proximity effect rely on the es- 
sential concepts in disordered Fermi-liquid physics, such 
as -Exh governing thermodynamic and quantum trans- 
port phenomena in mesoscopic systems. Since our Cm 
layer is a non-Fermi liquid, these concepts are not directly 
transferable. Therefore, to compare our findings with 
standard notions, we proceed along a phenomenological 
route frequently (but unwarrantedly) employed in exper- 
iments on SNS Josephson junctions with an unconven- 
tional N interlayer Namely, we extract an effective 
diffusion constant from the Kubo conductivity of the FK 
model |2^], using the Einstein relation upK = 2e^N{0)'D 
and the DOS at half-filling iV(0) from Fig. |, and then 
compute iJxh = KD/L^. It is interesting to check if such 
an energy scale provides any heuristic guidance in inter- 
preting our results at small t/pK where the hard mini- 
gap is present in the Cm spectrum. For example, in the 
largest Ufk = 2.0 sample, the resistivity of the FK cor- 
related metal is ppK — 240 /iilcm (assuming a — 3 A.). 
From here we get V « 2to^/fi. and £,n — 5.6a. This is sur- 
prisingly close to the true « 6.7a extracted from the 
decay of Ic in the SCmS junction as a function of the Cm 
layer thickness (the agreement improves for smaller 
Ufk)- The quasiclassical analysis for a mesoscopic dif- 
fusive junction of the same resistivity and with thickness 
L = 10a would give E'xh ~ O.IA and Eg « 0.32A. How- 
ever, at ?7fk = 2.0 no gap is found in the Cm spec- 
trum, while only a small dip (Fig. |^) in the LDOS per- 
sists as a remnant of the minigap opened for L/pK ^ 1-0. 
Moreover, the ratios of the minigap sizes E^/E^ in the 
SCmS junctions with two different thicknesses La — 5a 
and Lb = 10a are: (i) Ufk = 0.1 ^ E^/E''g = 2.9, 



(ii) Ufk 



0.25 



0.5 



E-'jEl = 3.1, (iii) Ufk 
A A, is a function of Ufk- These are different 
from the expected results of Eg/E^g = (Lb/La)^ = 4 or 
E^/E^g ^ Lb /La = 2, which would follow from Eg - E'xh 
analogy with the quasiclassical description of the prox- 
imity effect in a diffusive Fermi-liquid metal of the same 
resistivity [slJlll] , or a clean but chaotic interlayer , re- 
spectively. 

The last two pieces of information needed to charac- 
terize proximity-induced effects in correlated metals are: 

(1) the LDOS is position dependent, while Eg is spatially 
constant (which is the same as the quasiclassical phe- 
nomenology |ll|), and does not change upon lowering the 
temperature below our reference T/Tc = 0.091 for Fig. |; 

(2) the minigap is open for small enough [/pK, where An- 



dreev bound states also clearly coexist with it, as shown 
by the peaks below A in the Ufk = 0.1 case in Fig. |^. 
In fact, we find the largest minigap in the limit Ufk —> 
(corresponding to a clean SNS junction), which is gener- 
ated by the normal reflection at the SN interface due to 
a non-negligible A//i « 0.03 determining the amplitude 
of the scattering (/i is the Fermi energy measured 
from the bottom of the band). It appears that increasing 
of electron correlations by increasing Ufk then just leads 
to a monotonic vanishing of any initially open minigap in 
the noninteracting case [Eg{UFK = 0) depends on Ut, as 
illustrated in the inset of Fig. |l|, and is outside of the qua- 
siclassical approximation |2J] , but belongs to the realm of 
"noninteracting quasiparticle" physics]. This should be 
contrasted with the quasiclassical minigap as a function 
of quenched disorder [pl : the minigap forms for arbitrar- 
ily small concentration of impurities, increases with 1/i 
to a maximum value when i L, and then decreases in 
the diffusive limit as iJxh vpi- On the superconducting 
side of our SCmS junction we find that position depen- 
dent LDOS is nonzero in the energy range 2 A (Fig. ||), 
and decays to zero on the length scale of a few ^5 from 
the SCm boundary ||]. 
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